Wigner crystal and bubble phases in graphene in quantum Hall regime 
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Graphene, a single free-standing sheet of graphite with honeycomb lattice structure, is a semimetal 
with carriers that have linear dispersion. A consequence of this dispersion is the absence of Wigner 
[ crystallization in graphene, since the kinetic and potential energies both scale identically with density 

^v^j . of carriers. We study the ground state of graphene in the presence of strong magnetic field focusing 

on states with broken translational symmetry. Our mean-field calculations show that at integer 
fillings a uniform state is preferred whereas at non-integer fillings, Wigner crystal states (with broken 
translational symmetry) have lower energy. We obtain the phase diagram of the system. We find 
that it is qualitatively similar to that of quantum Hall systems in semiconductor heterostructures. 
' Our analysis predicts that non-uniform states, including Wigner crystal state, will occur in graphene 

in the presence of a magnetic field and will lead to anisotropic transport in high Landau levels. 

PACS numbers: 73.20.Qt, 73.43.-f 
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1 I. INTRODUCTION 

Two dimensional electron gas (2DEG), realized either in semiconductor heterostructures or by sprinkling electrons 
on Helium surface, has been extensively studied over the past few decades. It is known that such a system undergoes 
phase transition from a uniform state to a state with spontaneously broken translational symmetry, a Wigner crystal, 
when the carrier density is sufficiently low 1 or the temperature is sufficiently high. 2 ' In case of a 2DEG in semicon- 
' ductors, this state occurs when the gain from lowering the potential energy by localization outweighs the kinetic 
energy cost associated with the localization provided that the dispersion of carriers is well-described by an effective 
mass m* , Ek — ti 2 k 2 /2m*. This transition has been theoretically investigated^ although unequivocal experimental 
evidence is still lacking. The spectrum of the electron gas changes radically in the presence of a strong magnetic 
field. The kinetic energy is quantized and each Landau level (the manifold of eigenstates with a given energy) has 
| a macroscopic degeneracy. Since the kinetic energy is fixed for a given Landau level, it is possible to vary the ratio 
. of potential energy and kinetic energy by changing the filling factor within a given Landau level, and the system 
' undergoes a transition from a uniform state to states with broken translational symmetry. Indeed, this transition has 
been extensively studied theoretically^ There is strong experimental evidence that the ground state of such a system 
at filling factors v < 1/5 is a Wigner crystal^ and that at partial filling factors in high Landau levels, the ground 
state is non-uniform^ We remind the Reader that this evidence is obtained from transport measurements and that a 
direct measurement of spatial density modulation - crystalline structure - is exceedingly difficult since the 2DEG is 
buried under a substrate. 

In this paper, we focus on how these results change when the carries in the 2D gas have a linear dispersion instead 
of the the usual effective-mass quadratic dispersion. Graphene, a single sheet of graphite with honeycomb lattice 
structure, is a realization of a system with such carriers. It has the added advantage that such a 2D gas of carriers is 
not buried under a substrate and is, therefore, amenable to local probes that can investigate the crystalline structure. 
Graphene is a semimetal in which the valance and conduction bands touch at two inequivalent points K and K' = — K 
(and four other points related by symmetry). In the vicinity of these points (valleys), the band structure of carriers is 
well described by — ±hvak where Hvg = 5.8 eVA is the characteristic velocity and k is measured from one of the 
six points at which the conduction and valance bands touchji^ Due to the linear dispersion of carriers in graphene, 
potential and kinetic energies both scale as n 3 / 2 with the carrier density n. Therefore, graphene does not undergo 
Wigner crystallization in the absence of an external magnetic fields 

In the presence of a magnetic field, however, the kinetic energy is quantized and its ratio with the potential energy 
can be varied by changing the filling factor. This raises the prospect of Wigner crystal states in graphene. Here, 
we present a systematic mean-field analysis of the ground state of graphene in the quantum Hall regime, focusing 
on partial filling of the first few Landau levels. The plan for the paper is as follows. In the next section, we set up 
the low-energy Hamiltonian for graphene in the presence of magnetic field and recall results for the single-particle 
spectrum. In Sec. IIII1 we set up the formalism for Hartree-Fock (HF) approximation and discuss its details. In 
Sec. IIV1 we present the numerical results we obtain. We discuss the phase diagram of graphene as a function of the 
partial filling factor and the density profiles. We find that the results closely follow those of non-uniform states in the 
conventional 2DEG. We end the paper with a brief discussion and conclusions in section IVl 
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II. GRAPHENE IN MAGNETIC FIELD 

The low-energy Hamiltonian for electrons in the K valley is given b y 10 ' 11 

H K = V G (p X T X + PyTy) , (1) 

where t x and r y are Pauli matrices in the space consisting of two lattice sites A and B within a single unit cell 
(The Hamiltonian for the other valley is obtained by complex conjugation). The Hamiltonian in the presence of an 
magnetic field is obtained by Peirels substitution p > p — eA/c. In a uniform magnetic field B = Bz, generated by 
a vector potential A = Bxy, the Hamiltonian fl} becomes 

where Ck = — i[lBd x + (x/lB — fc?s)]/\/2 is the lowering operator, lg = ^JficjeB is the magnetic length, and [cfc,cjj = 1. 
The eigenvalues of this Hamiltonian are given by E n = ±hvay^\n\/lB and the corresponding eigenfunctions are given 
by 



(r|K,nfc) = -^= e lkv 



sgn(n)<£| n |_i(a; - kl%) 
<P\ n \(x - kl%) 



for n / and 





(r|K,0fc) = — l —e lky 



ip (x - kl B ) 



(3) 



(4) 



for n — 0. Here, ip n {x) are the simple harmonic oscillator eigenfunctions defined by <p n (x) — (x\{c\) n \0) / \fn\ and L y is 
the sample length in y-direction. (In the following, we will use units such that Is = 1)- Thus, n = eigenfunctions in 
graphene are identical to n = states in a conventional 2DEG, whereas for n^O, the eigenfunctions of graphene are 
an admixture of wavefunctions on the A and B lattice sitesJ^ Therefore, we expect that graphene, like conventional 
2DEG at partial filling factors, will support non-uniform (Wigner crystal) states. 

In the following, we denote crystals with one electron per unit cell, N e = 1 as Wigner crystals, and those with 
N e > 2 per unit cell as bubble crystals . 13 ' 14 We also consider modulated stripe states that can be described by oblique, 
rectangular, or centered rectangular latticesi 15 ' 16 We call these states anisotropic Wigner crystals, and the ones having 
a triangular or square lattice structure as isotropic Wigner crystals. 

III. HARTREE-FOCK HAMILTONIAN 

The microscopic Hamiltonian for carriers in graphene consists of the kinetic energy and Coulomb repulsion. We 
use a pseudospin notation to denote the valley index: a = + corresponds to the K valley and a = — corresponds to 
the K' = K valley. In the single-particle basis (|3I4[) the Hamiltonian is 

H = N^{E n ~ M°(0) + i Y, ^(^.^(^.^(-^^(-cOpS'.SW, (5) 

ncr B q,{ern} 

where iV^ = A/(2nl B ) is the number of flux quanta in the area A of the sample, fi is the chemical potential, 

V(q) = 2ne 2 /eq is the Coulomb interaction and e is the dielectric constant of graphene (e ~ 3), (q) is the density 
matrix element defined as 

-C''(q) = tJ- E e "* a " (k+fc ' ) ' asc L*wv«M'+* (6) 
^ kk' 

with c ank (c an k) the creation (annihilation) operator of the electrons, J 7 n n'{ < i} is the form factor 

(q) 

+ -e(\ ni \)e(\n 2 \) [F| ni | ) | na |(q)+Bgn(nin 2 )F| ni |_ 1) | n!1 |_ 1 (q)] , (7) 
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where 9{x) is the heavyside function and sgn(x) returns the sign of its argument. This form factor is a linear 
combination of form factors for wave functions on the two inequivalent lattice sites, 



F, 



n\ >n 2 



(q) 



1*2 1 

Iml 



{-% + iQx) 



V2 



(|ni| — |ri2 1) 



L 



|ni|-|n 2 | [ p -q 2 /4 
\n 2 \ 



(8) 



and F ni < n2 (q) = F n2 , ni (— q)*. Here L™\x) is the generalized Laguerre polynomial. To obtain the Hamiltonian |5|), 
we have ignored the interaction terms that scatter electrons from one valley to another and are exponentially and 
algebraically small in cl/Ib where a is the lattice constant of graphene (o ~ 5 A, Ib ~100 A)Ji 

The derivation of the HF Hamiltonian from Eq. ([5]) is straightforward and has been discussed extensively in the 
literature. 18 ' 19 ! 20 When inter-Landau level transitions are ignored, the Hartree-Fock Hamiltonian for a single Landau 
level index is given by 



Hhf — 



iy 2 



E 

ffQ 



{En ~ H) 

e 2 /d B 



p^(Q)-X^(Q)p^(Q) 



(9) 



where Q is a reciprocal lattice vector of the Wigner crystal and a = —a. The dimensionless Hartree and Fock 
potentials are given by 



=-Q 2 /2 



Q 



|^„,„(Q)| 2 pn(-Q)(l-'5Q,o), 



^r'(Q)= / dxe- x2 / 2 \F n , n (x)\ 2 J (xQ) P f(-Q) 
Jo 



(10) 
(11) 



where p^ a (Q) = (p^Zi (Q)) are determined self-consistently from the Hamiltonian (JOj) and p n (Q) = J2 a =± Pn' a (Q) 
is the total density at wavevector Q. 

The density matrix p^ a (Q) is determined from the equal-time limit (r — * _ ) of the single-particle Green's function 



Gr(fci,fc 2 ;r) = -<Tc nfclff (r)4 2(T ,(0)). 
We define the Fourier transform of G^ a (ki, fc 2 ; r) as 

G%°'(Q,iu m ) = ^ E f dre'^ k ^' 2+ ^ 6 k2 ^ Qy G° n °' (k u k 2 ;r) 
k lk2 Jo 

where (3 = 1/kgT is the inverse temperature. The equation of motion for G^ ,cr (Q, ioJ m ) is given b y 19 ' 20 
(iu m + p) 

where the Hartree and exchange self-energy terms are given by 







G„' + (Q,iw m ) 



E 

Q' 



£+>+(Q,Q') E+--(Q,Q') 
£-'+(Q,Q') S-.-(Q,Q') 



G+<+(Q',iw m ) 
G--+(Q',iw ro ) 



S^ CT (Q, Q') = [H n (Qf - Q) (Q' - Q)] exp [ -Q x Q'-l 



r (Q,Q') 



-X^(Q'-Q)exp( -QxQ'-l 



and Q x Q' ■ z = QxQ' v — QyQ' x - ln order to solve Eq. (fLT| , we diagonalize the self-energy matrix 



E 

Q' 



E+.+ (Q,Q') £+' 


-(Q,Q')" 




' W) " 




' Vj(Q) ' 


S-.+ (Q,Q') X-> 


-(Q,Q ; ). 






= w j 





(12) 
(13) 

(14) 

(15) 
(16) 

(17) 



in a basis with a specific lattice structure, where (Vj, Uj) is the j-th eigenvector and ujj is its corresponding eigenvalue. 
Using these eigenvectors and eigenvalues, one can calculate the density matrix 



Pn' + (Q) Pn'~(Q) 



J2[V*(0),U* 



Uj(Q) 



f(uj ~ P) 



(18) 
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where f(u> — fj,) is the Fermi-Dirac distribution function and the chemical potential fi is determined by 

P^ + (0) = E Vj(0)V;(0)f(^ - M) = \v (19) 

3 

where v is the partial filling factor in the Landau level n. We solve the set of equations (| 1 441 1 9 [) self-consistently 
to obtain the resultant density matrix and the ground state energy per particle for different lattice structures The 
real-space density profile is then obtained using inverse Fourier transform, 

= ^pE r(Q)^,„(Q)e lC! ' r (20) 

B crQ 

For simplicity, in the numerical results for density profile, we use the dimensionless density, 2-7r^n(r). 



IV. NUMERICAL RESULTS 

Now we turn to the numerical results. It follows from Eq. ([9]) that a state with intervalley coherence always has 
lower energy than a state without the coherence. Therefore, we only focus on solutions with Pn'~(Q) 0- Thus, 
the ground state consists of electrons occupying the symmetric state (between the two valleys) when the filling factor 
within a Landau level is v < 1. (Here, we have assumed spin-split Landau levels- 1 and ignored inter-Landau level 
transitions because they do not qualitatively change our conclusions.) 

We use a simplified oblique lattice with two primitive lattice vectors^ ai = (a, 6/2), a.2 = (0, b) (We perform similar 
analysis with square, rectangular, and centered rectangular lattices, and obtain results similar to those presented 
below). We denote the ratio 7 = b/a. Note that the triangular lattice (7 = 2/y/S) and quasi-striped states (7 — > 0) 
are its special cases. We do not consider purely one-dimensional striped states because they are prone to density 
modulations along the stripes^ The tw o lattice c onstants are determined by the constraint that the unit cell contains 
N e electrons, and are given by a — ls\/2TrN e /^ and b = aj. The reciprocal lattice basis vectors are 

2vr , N , 2tt / 1 1 \ 
bx = — (1,0), b 2 = — --,- , (21) 



a a 



2' 7 



and the reciprocal lattice vectors are given by Q mn = mh\ + nlo2- We determine the optimal lattice structure by 
obtaining the 7 (0 < 7 < 2/y/S) that minimizes the mean-field energy. We use an external (infinitesimal) potential 
with the symmetry of the lattice to generate the initial density matrix, 



£ e x(Q, Q') = -exp 



-f^(Q 2 + Q' 2 ) + ^(Q-Q' + *QxQ'-z) 



(22) 



where Qq = |Qn|. We adjust the number of basis states used to calculate the density matrix, and verify that the 
zero-temperature sum rulei2i 

E [l^ + (Q)| 2 + \Pn'-(Q)\ 2 ] = P+' + (0) (23) 
Q 

is satisfied within an accuracy of 10~ 5 . 

For n — 0, the equivalence between single-particle wavefunctions for graphene and conventional 2DEG, implies 
that at small filling factors, the ground state is a triangular Wigner crystal^ Indeed, graphene, with its pseudospin 
(valley index) can be mapped onto a bilayer system in which layer index is the pseudospin, in the limit when the 
layer separation d — > 0. Our calculations reproduce the results for ground state energy and lattice structure. Figure 
[T] shows the mean-field energy per particle (phase diagram) as a function of partial filling factor for different lattice 
structures. We see that for v < 0.62, the ground state is triangular Wigner crystal; it becomes an anisotropic Wigner 
crystal for v > 0.62. We find that a bubble crystal with N e = 3 is identical, in energy, to the triangular Wigner 
crystal, and the bubble crystals with N e > 2 have higher energies. Figure [5] shows the real-space electron density 
profile for graphene at partial filling v = 0.25 and for a bilayer quantum Hall system at layer separation d = at the 
same filling factor. We obtain, as expected, identical density profiles with a triangular Wigner crystal. As the partial 
filling v is increased, the ground state of the system changes to an anisotropic Wigner crystal (Figure [3]). We find, in 
general, that the optimal value of anisotropy is high. Therefore, the electron density resembles uniform stripe states 
and the density modulation along the stripes is quite small (Figure [3]). 
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FIG. 1: (Color Online) Ground state energy per particle (measuyed in units of e 2 /els) in graphene as a function of partial 
filling v in the n — Landau level. The ground state is a triangular Wigner crystal for small v and becomes an anisotropic 
Wigner crystal as v increases. The arrow mark critical filling factors v* at which transition from the triangular Wigner crystal 
to the anisotropic Wigner crystal takes place. Note that for n — 0, the bubble crystal always has higher energy. 




FIG. 2: (Color Online) A density plot of the ground state electron density in graphene (left) and a bilayer system with d — 
(right) in the n — Landau level with partial filling v = 0.25. According to the phase diagram in Fig. [1] the ground state is 
an isotropic triangular crystal. As expected, due to the equivalence between the two systems in the n = Landau level, the 
two densities are identical. 

For n > 1, the equivalence between graphene and a bilayer quantum Hall system at layer separation d = breaks 
down since the single-particle wave functions are different. Figures |4][6] show the phase diagram as a function of partial 
filling factor v for Landau levels n=l-3, respectively. The phase diagram of n = 1 Landau level is qualitatively similar 
to n — Landau level. For n > 2 we find that the ground state for graphene is an isotropic crystal when the partial 
filling factor v is small, whereas it is an anisotropic crystal when v is sufficiently large. For intermediate values of v, 
we find that the ground state is a bubble crystal with N e = 2 or N e — 3. For example, for n = 2 Landau level in 
graphene, the triangular Wigner crystal is stable for v < 0.28, the bubble state with two electrons N e = 2 is stable 
for 0.28 < v < 0.43, whereas for v > 0.43 an anisotropic Wigner crystal has the lowest energy. We can also see 
that the critical values of v at which transitions from a triangular Wigner crystal to a bubble state to an anisotropic 
Wigner crystal take place are systematically higher than corresponding values for a bilayer system at d = 0. Figure 
[7] shows optimal value of 7(f) for the ground state crystal structure for different Landau level indices. We see that 
the transition from an isotropic Wigner crystal (7 = 2/V3 = 1.15) to an anisotropic Wigner crystal in graphene (a) 
occurs at higher values of v than the corresponding values in bilayer systems (b). 

Figure [8] shows the (dimensionless) real-space electron density profile for the n = 3 Landau level when the system 
is in the isotropic Wigner crystal state (y = 0.18), in the bubble crystal state with N e = 2 (y = 0.25) and with 
N e = 3 (y = 0.35), and in the anisotropic Wigner crystal state [y = 0.75). We see that the density profile is different 
from that in the lowest Landau level, due to different form factors. It is instructive to compare density profiles of the 
anisotropic Wigner crystal in the n — 3 Landau level and n = Landau level (Fig. [3|). We see that both resemble a 
(quasi)-uniform striped states, although the n — 3 density profile also shows the existence of uniform stripes between 
the modulated ones. 

We end this section by comparing the electron density profiles in graphene and bilayer system with d — at partial 
filling v = 0.25 in the n = 1 Landau level in Fig. [9l Recall that in the lowest Landau level, these density profiles are 
identical (Fig|2J). For both systems, the self-consistent solutions for the density matrix p^ 7 (q) are identical; however, 
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FIG. 3: (Color Online). Electron density profile for the anisotropic Wigner crystal ground state at v — 0.75 in the lowest 
Landau level. Notice that the density modulation along the y-axis is small and the system behaves essentially like a striped 
state. 




FIG. 4: (Color Online) Ground state energy per particle (measured in units of e 2 /eZs) for different crystal structures in the 
n = 1 Landau level for (a) graphene, and (b) a bilayer system at d = 0. The arrows mark the critical filling factor u* at which 
transition from the triangular Wigner crystal to anisotropic Wigner crystal takes place. 



for n > 1, due to the differences in form factors, the resultant real-space electron density profiles are different. 



V. DISCUSSION AND CONCLUSIONS 



In this paper, we have systematically studied the ground state of graphene in the presence of a strong magnetic 
field, focusing on broken symmetry states with intervalley coherence, using the Hartree-Fock mean-field analysis. We 
have ignored the inter-Landau level transitions (since their inclusion does not qualitatively change the phase diagram) 
and considered spin-split Landau levels. We have focused only states with inter- valley coherence because, due to the 
isotropic Coulomb interaction, energy of a state without coherence is always larger than the energy of a state with 
intervalley coherence. Thus, in experiments, these levels will correspond to total filling factors < vt < 1 (lowest 
subband of the n = level), 2 < vt < 3 (lowest subband of the n = 1 level), 6 < vt < 7 (lowest subband of the n = 2 
level). Our calculations, naturally, have also not taken into account competing (Laughlin) fluid states with uniform 
density. Such fluid states will have lower energy at special filling factors, for example, v = 1/3; however, our results 
likely represent the true ground state of the system at generic filling factors and in high Landau levels (where, in a 
conventional 2DEG, Hartree-Fock solutions are reliable). 
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FIG. 5: (Color Online) Ground state energy per particle (measured in units of e 2 /e/s) for different crystal structures in the 
n = 2 Landau level for (a) graphene, and (b) a bilayer system at d — 0. The arrows mark critical filling factors v* at which 
transition from the triangular Wigner crystal to bubble crystal to anisotropic Wigner crystal take place. 




FIG. 6: (Color Online) Ground state energy per particle (measured in units of e 2 /e/s) for different crystal structures in the 
n = 3 Landau level for (a) graphene, and (b) a bilayer system at d — 0. The arrows mark critical filling factors v* at which 
transition from the triangular Wigner crystal to bubble crystal to anisotropic Wigner crystal take place. 



Our analysis found that, since the kinetic energy of graphene's linearly dispersing carriers is quenched in the 
presence of a magnetic field, graphene is qualitatively similar to a conventional 2DEG with a quantum degree of 
freedom (bilayer quantum Hall system at d = or a single layer system with vanishing Zeeman coupling) . We showed 
that the (mean-field) ground state of graphene is evolves from a triangular Wigner crystal to an anisotropic Wigner 
crystal (striped state) as the partial filling v in a given Landau level is increased. We also showed that for Landau 
level indices n > 2, a bubble crystal state with two electrons per unit cell, N e = 2 occurs at intermediate values of 
We have compared our results with mean-field results for a corresponding bilayer system at d = (where the 
Coulomb interaction becomes isotropic in the pseudospin space). Our findings indicate that different form factors 
for electrons in graphene systematically shift the critical values of v at which the phase transitions occur to higher 
values, thus expanding the region of stability for the triangular Wigner crystal and bubble crystal states, compared to 
their counterparts in bilayer systems. 

These broken symmetry states in graphene are open to more probes than the conventional 2DEG. Our analysis 
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FIG. 7: (Color Online) Optimal crystal structure parameter "y(v) for the ground state of (a) graphene and (b) bilayer system 
at d — (bottom) for different Landau level indices. Recall that 7 = 2/\/3 = 1.15 corresponds to a triangular lattice whereas 
7 < 0.5 corresponds to highly anisotropic Wigner crystals or striped states. 




FIG. 8: (Color Online) A density plot for the ground state electron density in graphene with Landau level n = 3. Shown here 
are a triangular Wigner crystal at v = 0.18 (top left), a bubble state with N e = 2 at v — 0.25 (top right), a bubble state with 
N e — 3 that occurs at v = 0.35 (bottom left), and the anisotropic Wigner crystal state at v = 0.75 (bottom right). 

predicts that anisotropics in the longitudinal resistance will be observed in high Landau levels (similar to those 
observed in quantum Hall systems) and will provide a signature of highly anisotropic Wigner crystal (striped) states. 
In addition, since the 2DEG in graphene is literally at the surface, the electron density modulations can be directly 
probed, by scanning tunneling microscopy, and may provide a direct evidence of electronic crystal states. A weak 
random disorder will, in general, pin these crystals and destroy the true long-range order. However, since it does 
not prefer one crystal structure over another, the disorder will not change the phase diagram qualitatively. Direct 
experimental observation of charge and current density distribution in quantum Hall systems is an outstanding 
problem; observation of isotropic Wigner crystal and striped states in graphene will improve our understanding of 
local structure of quantum Hall states. 
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FIG. 9: (Color Online) A density plot of the ground state electron density profile in graphene (left) and bilayer system at 
d — (right) at partial filling factor v = 0.25 in the n = 1 Landau level. Both systems have identical density matrices p n (Q). 
However, the difference in the form factors in Eq. (|20p leads to different real space densities. 
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